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Abstract 



^ • We summarize the published top-quark mass measurements from the CDF and D0 ex- 

^ periments at Fermilab. We combine published Run-I (1992-1996) measurements with the 

most recent published Run-II (2001-present) measurements using up to 340 pb~ of data. 
Taking correlated uncertainties properly into account the resulting mass of the top quark 
is Mt = 174.2 lb 2.0(stat) ± 2.6(syst) GeV/c^, which corresponds to a total uncertainty of 
3.3 GeV/c^, i.e. 1.9% precision. 

Since this combination uses only a subset of the available analyses and data sets, it 
does not supersede our latest world average combination of Mt = 172.5 ± 1.3(stat) ± 
1.9(syst) GeV/c^, which is based on the latest published and preliminary results. 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton colhder located 
at the Fermi National Accelerator Laboratory, have made several direct experimental measure- 
ments of the top-quark pole mass, Mt. The pioneering measurements were based on about 
100 pb'^ of Run-I (1992-1996) data El IHjSl El El El El CHI CIl [EI and include results 
from the tt — >• qq'bqq'b (all-j), the tt — > iuqq'bb (}+]), and the ti -^ i~^i'bi~vb (di-1) decay 
channels^. Results using approximately 340 pb~ of Run-II (2001-present) data have been re- 
cently published in the l-|-j and di-1 channels P^ ITU [T3] . More preliminary analyses have been 
performed in the l-|-j and di-1 channels using 370 — 750 pb^^ of data and improved analysis 
techniques |inilI7l[IHllinil2niEIll22ll2Sll21l2nil2Sll2Z|- 

This note reports the average top quark mass obtained by combining five published Run-I 
measurements [21 El El El QHl Q] with the most recent published Run-II measurements from 
CDF ^3 El El -^ Since this combination only uses a subset of the available analyses and data 
sets, it does not supersede our latest world average combination of Mt = 172.5 ± 1.3(stat) ± 
1.9(syst) GeV/c^ J2H1- The combination takes into account the statistical and systematic un- 
certainties and their correlations using the method of references |29| EDj- The most precise 
individual measurements of M^ in this combination are the measurements in the l-|-j channel 
from Run I and Run 11. These are 173.5l|l GeV/c^ (CDF II, ^) and 180.1 ± 5.3 GeV/c^ 
(D0 I, [ID])- These have weights in the new Mt combination of 57% and 26%, respectively. 

The input measurements and error categories used in the combination are detailed in Sec- 
tion El and El respectively. The correlations used in the combination are discussed in Section JH 
and the resulting world average top-quark mass is given in Section El A summary and outlook 
are presented in Section El 



2 Input Measurements 

For this combination seven measurements of Mt are used, five published Run-I results and two 
published Run-II results. The most precise result in each channel is considered for each experi- 
ment. In general, the Run-I measurements all have relatively large statistical uncertainties and 
their systematic uncertainty is dominated by the total jet energy scale (JES) uncertainty. In 
Run-II both CDF and D0 take advantage of the larger ti samples available and employ new 
analysis techniques to reduce both these uncertainties. In particular the JES is constrained 
using an in-situ calibration based on the invariant mass oi W ^ qq' decays in the 1+j channel. 



^Hcre £ = e or /i. Decay channels with exphcit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 

•^D0 results in the l+j and di-1 channels using Run II data are at present close to submission but not yet 
published and are therefore not included in this combination. 



The Run-II CDF analysis in the 1+j channel constrains the JES using the in-situ W -^ qq' 
decays. The JES is also determined using "external" calibration samples as was done for the 
Run-I measurements. This external JES is applied as an additional constraint in this CDF 
analysis to improve further the total JES uncertainty. Small residual JES uncertainties arising 
from rj— and p^-dependencies and the modeling of 6-jets are included in separate error cate- 
gories. The Run-II CDF measurement in the di-1 channel uses only the externally determined 
JES, some parts of which are correlated with the externally determined Run-I JES as noted 
below. 

The CDF Run-II measurement in the l-|-j channel requires special treatment in order to ac- 
count more accurately for the correlations of the JES uncertainties since the fit uses information 
from both the in-situ and external JES calibrations. In the combination we treat this one mea- 
surement as two separate inputs - one which includes only the in-situ JES calibration, (l+j)i, 
and a second which includes only the JES as determined from the external calibration samples, 
(l-|-j)e. We correlate the JES related error categories as described below while taking the rest of 
the error categories to be 100% correlated between these two inputs. The combination of just 
these two inputs using these correlations yields the identical central value, statistical, JES, and 
total systematic uncertainty as the measurement reported in reference ^H]- The correlations 
between these two inputs and the rest of the inputs are as described in Section HI 

The inputs used in the combination are summarized in Table C] with their uncertainties 
sub-divided into the categories described in the next Section. 



3 Error Categories 

We employ the same error categories as used for the previous world average [31]. They have 
evolved to include a detailed breakdown of the various sources of uncertainty and aim to lump 
together sources of systematic uncertainty that share the same or similar origin. For example, 
the "Signal" category discussed below includes the uncertainties from ISR, FSR, and PDF - 
all of which affect the modeling of the tt signal. Additional categories have been added in 
order to accommodate specific types of correlations. For example, the jet energy scale (JES) 
uncertainty is sub-divided into several components in order to more accurately accommodate 
our best estimate of the relevant correlations. Each error category is discussed below. 



Statistical: The statistical uncertainty associated with the Mt determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 





Run-I published 


Run-II published 




CDF 




D0 




CDF 






all-j 


1+J 


di-1 


1+J 


di-1 


(1+J). 


(l+j)e 


di-1 


Lumi (pb ^) 


110 


105 


110 


125 


125 


320 


320 


340 


Result 


186.0 


176.1 


167.4 


180.1 


168.4 


173.5 


173.5 


165.2 


iJES 


0.0 


0.0 


0.0 


0.0 


0.0 


4.2 


0.0 


0.0 


bJES 


0.6 


0.6 


0.8 


0.7 


0.7 


0.6 


0.6 


0.5 


cJES 


3.0 


2.7 


2.6 


2.0 


2.0 


0.0 


2.0 


2.2 


dJES 


0.3 


0.7 


0.6 


0.0 


0.0 


0.0 


0.6 


0.8 


rJES 


4.0 


3.4 


2.7 


2.5 


1.1 


0.0 


2.2 


1.1 


Signal 


1.8 


2.6 


2.8 


1.1 


1.8 


0.8 


0.8 


1.3 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.5 


0.5 


0.8 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.6 


0.6 


1.3 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.2 


0.2 


0.8 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.8 


3.9 


3.6 


4.4 


3.3 


3.4 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


2.7 


2.7 


6.1 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


5.2 


4.3 


7.0 



Table 1: Summary of the measurements used to determine this average M^ using only published 
Mt results. As described in the text, the CDF Run-II measurement in the lepton+jets channel 
is treated as two inputs in order to more accurately account for the correlations in the jet energy 
scale uncertainties. All numbers are in GeV/c^. The error categories and their correlations are 
described in the text. The total systematic uncertainty and the total uncertainty are obtained 
by adding the relevant contributions in quadrature. 



W — > qq' invariant mass in the CDF Run-II 1+j measurement. Residual JES uncertain- 
ties, which arise from effects not considered in the in-situ calibration, are included in 
other categories. 

bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modeling of 6-jets and which is correlated across all measurements. For both CDF and 
D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, 6-fragmentation modeling, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run-II studies but back-propagated to the 
Run-I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 



cJES: That part of the JES uncertainty which originates from modehng uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibra- 
tion data samples used and which is correlated between measurements within the same 
data-taking period (i.e. Run I or Run II) but not between experiments. For CDF this 
corresponds to uncertainties associated with the //-dependent JES corrections which are 
estimated using di-jet data events. 

r JES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. This is dominated by uncertainties in the calorimeter response to light- 
quark jets. For CDF this also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the ti signal 
which is correlated across all measurements. This includes uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the tt Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modeling the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 
uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
l-|-j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W+jets 
background (all channels) are included. 

Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to calibrate each method. 

Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
includes variations observed when substituting PYTHIA [^ EHl E3 (Run I and Run II) 
or ISAJET [SS] (Run I) for HERWIG jHEl EZ] when modehng the ti signal. Similar 
variations are included for the D0 Run-I measurements. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in 
the D0 calorimeter and from the multiple interaction corrections to the JES. 



These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 



method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run-I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (<^ 100 MeV/c^) the combined M^ and total uncertainty. 



4 Correlations 



The following correlations are used when making the combination: 

• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrected 
among the measurements. 

• The uncertainties in the dJES category is taken to be 100% correlated among all Run-I 
and all Run-II measurements on the same experiment, but uncorrelated between Run I 
and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 

• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 

Using the inputs from Table Q] and the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Table |21 

The measurements are combined using a program implementing a numerical x^ minimiza- 
tion as well as the analytic BLUE method j^niEOl- The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination jHS], and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments puj . 



5 Results 

The combined value for the top-quark mass is: 

Mt = 174.2 ± 3.3 GeV/c^ , (1) 
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Run-I published 


Run-II published 




CDF 




D0 


CDF 






1+J 


di-1 


all-j 


1+J 


di-1 


(1+j). (l+j)e di-1 


CDF-I 


1+J 


1.00 












CDF-I 


di-1 


0.29 


1.00 










CDF-I 


all-j 


0.32 


0.19 


1.00 








D0-I 


1+J 


0.26 


0.15 


0.14 


1.00 






D0-I 


di-1 


0.11 


0.08 


0.07 


0.16 


1.00 




CDF-II 


(1+J). 


0.08 


0.05 


0.03 


0.07 


0.03 


1.00 


CDF-II 


(l+j)e 


0.51 


0.29 


0.34 


0.26 


0.11 


0.41 1.00 


CDF-II 


di-1 


0.26 


0.17 


0.18 


0.17 


0.09 


0.04 0.30 1.00 



Table 2: The resulting matrix of total correlation coefficients used in the top-quark mass 
combination reported here. 

with a x^ of 5.8 for 6 degrees of freedom, which corresponds to a probability of 45% indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
into the contributions from the various error categories as: Statistical (±2.0), total JES (±2.3), 
Signal (±1.0), Background (±0.6), Fit (±0.4), Monte Carlo (±0.2), and UN/MI (±0.4), for a 
total Systematic (±2.6), where all numbers are in units of GeV/c^. The pull and weight for 
each of the inputs are listed in Table 121 The input measurements and the resulting combined 
to-quark mass are summarized in Figure ^ 



In this combination, for the CDF Run-II l±j measurement using 320 pb~ of data, the in-situ 
and external JES calibrated inputs each carry approximately the same weight. In the latest 
world average combination reported Reference |2H]; the weight of the CDF Run-II 1+j input 
using the in-situ JES calibration carries three times the weight of its counterpart using the 
external JES calibration. This trend is expected to continue with more data since the in-situ 
JES uncertainty is expected to improve as the statistics of the W —>■ qq' sample increase with 
larger data sets. In contrast the uncertainty on the external JES calibration already has large 
contributions from modeling uncertainties which may not be reduced with larger data sets. 



Although the x^ from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, l±j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 



rall- 



rl+J 



fdi-l 



physical observables, M^ , M^ , and M^ ' . The results of this combination are shown in 
Table Eland have x^ oi 1.4 for 5 degrees of freedom, which corresponds to a probability of 93%. 





Run-I published 


Run-II published 




CDF 


D0 


CDF 




1+J 


di-1 all-j 


1+j di-1 


(1+j). (l+j)e di-1 


Pull 


+0.30 


-0.62 +1.07 


+ 1.43 -0.46 


-0.17 -0.24 -1.45 


Weight [%] 


1.0 


1.1 0.6 


26.2 2.8 


28.2 28.5 11.6 



Table 3: The pull and weight for each of the inputs used to determine this average mass of the 
top quark. 



Parameter 


Value (GeV/c2) 


Correlations 


^all-j 


186.7+11.0 
175.6+ 3.4 
166.1+ 5.7 


1.00 

0.27 1.00 

0.16 0.35 1.00 



Table 4: Summary of the combination of the seven measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton+jets, and di-lepton 
channel. 

These results differ from a naive combination, where only the measurements in a given channel 
contribute to the M^ determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. 



6 Summary 

A combination of published measurements of the mass of the top quark from the Tevatron exper- 
iments CDF and D0 is presented. The combination includes five published Run-I measurements 
and two published Run-II measurements. Taking into account the statistical and systematic un- 
certainties and their correlations, the average top quark mass result is: Mt = 174.2+3.3 GeV/c^. 



This average of published results has an accuracy of 1.9%, while the present world average 
top quark mass of Mt = 172.5 ± 2.3 GeV/c^ J2H], which also includes preliminary analyses 
of more data by CDF and D0, has an accuracy of 1.3%. Both combinations are limited by 
the systematic uncertainties, which in turn are dominated by the jet energy scale uncertainty. 
This systematic is expected to improve as larger data sets are collected since new analysis 
techniques constrain the jet energy scale using in-situ W —>■ qq' decays. It can be reasonably 
expected that with the full Run-II data set the top-quark mass could be known to better than 
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1%. To reach this level of precision further work is required to determine more accurately the 
various correlations present, and to understand more precisely the 6-jet modeling, Signal, and 
Background uncertainties which may limit the sensitivity at larger data sets. Limitations of 
the Monte Carlo generators used to calibrate each fit method may also become important as 
the precision reaches the 1% level and will warrant further study in the future. 



Mass of the Top Quark 



Measurement 
CDF-I di-l - 
D0-I di-l — 
CDF-II di-l 
CDF-I l+j 
D0-I l+j 
CDF-II l+j 
CDF-I all-j 



Tevatron Run- 



150 



M 



M^^p [GeV/c1 
167.4 ±11.4 
168.4 ±12.8 
165.2± 7.0 
176.1 ± 7.3 
180.1 ± 5.3 
173.5± 3.9 

- 186.0 ±11.5 



X /dot = 5.8/6 



174.2± 3.3 



tOD 



170 
[GeV/c^] 



190 



Figure 1: A summary of the input measurements and the resulting combined top-quark mass. 
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